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Abstract 

The Faddeev-Volkov solution of the star-triangle relation is connected with the modular 
double of the quantum group U q (sl2)- It defines an Ising-type lattice model with positive 
Boltzmann weights where the spin variables take continuous values on the real line. The free 
energy of the model is exactly calculated in the thermodynamic limit. The model describes 
quantum fluctuations of circle patterns and the associated discrete conformal transformations 
connected with the Thurston's discrete analogue of the Riemann mappings theorem. In 
particular, in the quasi-classical limit the model precisely describe the geometry of integrable 
circle patterns with prescribed intersection angles. 
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1 Introduction 



The Faddeev-Volkov solution [1-3] of the Yang-Baxter equation [4] possesses many remarkable 
properties. From algebraic point of view it is distinguished by its modular duality. The solution 
is related with non-compact representations of the modular double [5] of the quantum group 
U q {sl2) <8> Uq^sli), where q = e l7rft2 and q = e _i7r / fe2 . It was studied in connection with the lattice 
Liouville and sinh-Gorgon models [6-9]. Note, that the parameter b above is related to the 
Liouville central charge cl = 1 + 6 (b + b^ 1 ) 2 . 

It is well known that every solution of the Yang-Baxter equation defines an integrable model 
of statistical mechanics on a two-dimensional lattice. The Faddeev-Volkov model, defined by 
their solution, has strictly positive Boltzmann weights. It is an Ising-type model with continuous 
spin variables cij 6 R, taking arbitrary real values. Its partition function is given by the integral 

Z= I e" £ M Y[da h (1) 

J i 

where the energy £[p] is a convex and bounded from below function of the spin variables a = 
{o~i, o~2, ■ ■ •}■ In this paper we study this model and reveal its remarkable connections to discrete 
geometry, namely, to circle patterns and discrete conformal transformations (for an introduction 
see a popular review [10] and a monograph [11]). Recently these topics attracted much attention 
due to Thurston's circle packing approach to the discrete Riemann mapping theorem [12-15] 
and the discrete uniformization theorem [16-19]. 

In this approach circle packings and, more generally, circle patterns (see Sect. 5 below) serve 
as discrete analogues of conformal transformations, where the circle radii can be thought as the 
local dilatation factors. They satisfy certain integrable non-linear difference equations which 
play the role of discrete Cauchy-Riemann conditions. These equations are called the cross- 
ratio equations [20-22]. In the simplest case they reduce to the Hirota form of the discrete 
sine-Gordon equation [23,24]. (Curiously enough, the continuous Riemann mapping problem 
for domains with a piece- wise smooth boundary is also connected to the Hirota equations [25], 
but of the different type [26].) An appropriate variational approach for the circle patterns was 
developed by Bobenko and Springborn [27]. Their action A[p] is a positive definite function 
of the logarithmic radii of the circles p = {pi, pi . . .}, Pi = logrj. We show that the leading 
quasi-classical asymptotics of the energy functional in (pQ) 

E[a]-^-^A[2TTba), 6^0, c L -► +oo, (2) 

precisely coincide with this action, provided logarithmic radii of the circles are identified with 
the (rescaled) spin variables in ([1]), namely p% = 2irbai. Thus, at cl — ► +oo, stationary points of 
the integral (P) correspond to circle patterns. At finite values of cl the model describes quantum 
fluctuations of these patterns and the associated discrete conformal transformation. Given that 
the circle radii define local dilatation factors, the Faddeev-Volkov model describes a quantum 
discrete dilaton. 

The continuous quantum field theory with the conformal symmetry [28] has remarkable 
applications in physics and mathematics [29]. It would be interesting to understand which 
aspects of the continuous conformal field theory can be transferred to the discrete case. 

The organization of the paper is as follows. In Sect. 2 we present definitions and discuss 
properties of the Faddeev-Volkov model. In Sect. 3 we consider Z-invariant lattices and rhombic 
tilings which are closely related to the combinatorics of the circle patterns. The quasi-classical 
limit is considered in Sect. 4. An introduction to circle patterns is given in Sect. 5. We review the 
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derivation of the cross-ratio equations and show that they are exactly coincide with the quasi- 
classical equation of motion in the Faddeev-Volkov model. Some open problem are discussed in 
the Conclusion. In Appendix A we discuss some implications of the star-triangle equation and 
.Z-invariance to the circle patterns and volumes of polyhedra in the Lobachevskii 3-space. The 
Appendix B contains useful information on the special functions used in this paper. 

2 The Faddeev-Volkov model 

Consider an Ising-type model on the square lattice, shown in Fig(TJ Each site i of the lattice is 
assigned with a spin variable <7j G R, taking continuous real values. Two spins a and b interact 
only if they are adjacent (connected with an edge of the lattice). Typical horizontal and vertical 
edges are shown in Figj2j The corresponding Boltzmann weights are parametrized through the 




Figure 1: The square lattice (solid lines) and its medial lattice (dashed lines). 

so-called "rapidity variables" p and q associated with the oriented dashed lines in Fig(TJ In our 
case the weights depends only on the spin and rapidity differences a — b and p — q, where a and 
b are the spins at the ends of the edge. We will denote them as W p - g (a — b) and W p - q {a — b) 
for the horizontal and vertical edges, respectively. 
The partition function is defined as 

z = [ n w p-m - ^ n w p-i( a k - <ri) n ^ . (3) 

(m) (fc,0 

where the first product is over all horizontal edges the second is over all vertical edges 

(k,l). The integral is taken over all the interior spins; the boundary spins are kept fixed. 

Explicit expressions for the Boltzmann weights contain two special functions. The first one, 

I \ d * 1 l\ [ e ~ 2 ' ZW ^\ (A\ 

<p[z) - exp ^ sinHwb)sinh{w/b) w ) > W 
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Figure 2: Typical horizontal and vertical edges and the corresponding Boltzmann weights. 

is the Barnes double-sine function [30]. It was introduced to this field in [3] as the non-compact 
version of the quantum dilogarithm [31]. The second function has no special name, 

, , def (1 [ e~ 2lzw dw\ 

{Z) " 6XP VsyK+iosinh^sinh^-^cosh^^ + fe- 1 )) w) ' [ ) 

but previously appeared [32] in the context of the quantum sinh-Gordon model. The properties 
of these functions are briefly summarized in the Appendix IB1 
The Boltzmann weights W$ (s) and Wg (s) are defined as 

We (s) d ^ -L , We (s) * W^s) , (6) 

where 9 and s stand for the rapidity and spin differences, respectively, and the parameter 77 
reads 

r ] = (b + b~ 1 )/2. (7) 
The normalization factor Fe has the form 

Fe <M e i^ 2/ ^ (1 _ 8r? 2 )/24 $(2i7?0/7r) (8) 

The positivity requirement for the Boltzmann weights (point (b) below) singles out two different 
regimes 

(i) b > 0, or (ii) |6| = 1, Im(6 2 ) > 0. (9) 

The following consideration apply to either of these cases. Note that due to the symmetry 
b <-> for the regime (i) it is enough to consider the range < b < 1. The weights ([6]) possess 
the following key properties: 

(a) Symmetry: Wg(s) and We{s) are even functions of s. The definition © implies 

W e (s) = W„- e (s) . (10) 

(b) Positivity: when the parameter b belongs to either of the regimes ([9]) the functions Wg (s) 
and We{s) are real and positive for < 9 < it and real s . They have no zeros and poles on 
the real axis of s. 

(c) Asymptotics: 

W e {s) ~ F -V 2 ^l s l(l + O(e- 27r/3 l s l)), s -> ±00 . (11) 
where j3 = min(Re&, Reft -1 ). 
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(d) Star-triangle relation: 

/ daW q - r (a-a)Wp-r(c-a)W p - q (a-b) = W p - q (c-a)W p - r (a-b) W q - r (c-b) . (12) 

(e) Inversion relations: 

W e {a-b)W„ e {a-b) = 1, lim / dc W H+£ (a - c)W .- d+e (c - b) = 5(a-b), (13) 



where t is real, 
(f) Self-duality: 



W 9 (s) = I dx e 2nixs W 9 (x) . (14) 



(g) Initial conditions: 

W e (s) = 1 + 0(0) and W e {s) = 5{s) + 0{9) as 9^0, (15) 
where S(s) is the <5-function. 

The normalization factor Fe is the "minimal" solution of the inversion and crossing symmetry 
relations, [33-35] 

FdF ^ e = li J±_ = e -2^0(,-0)/,+ m (i+^yn (?W_ ir] \ (16) 
F-K-e V n ) 

in the sense that it does not have poles and zeroes in the strip — ir < He 6 < it; it is real and 
positive for real —it < 9 < it, and has a simple pole at 9 = tt and a simple zero at 9 = —tt. Note 
also that F = 1. Equation (fT4"j) is a consequence of (jB.8[) and the second relation in (fTUj) . 

Let AT be the total number of edges on the lattice. In the thermodynamic limit, when N is 
large, the partition function grows as 

log Z = -f bu lkN - fboundary^/N + 0(1), N -> OO, (17) 

where the first and second terms describe the bulk and boundary contributions to the free energy. 
The inversion relation arguments of [33-35] imply that with our normalization of Boltzmann 
weights ([6]) the bulk free energy for the system ([3]) vanishes exactly 

-f bu ik= lim -J-logZ = (18) 



Thus, log Z can only grow linearly in the number of boundary edges, which is 0(\/N). Here 
we assume that the boundary spins are kept finite in the limit N — > oo, otherwise (|18p will not 
hold. 

The weights (|6|) attain their maxima at s = 0. In the quasi-classical domain < b < 1 the 
constant W* (0) slowly varies between the values 

Wn(0) = e? = 1.3385. .. , Wn(0) = V2 = 1.4142 ... , (19) 

2 b=0 2 b=l 

where G = 0.915965 ... is the Catalan's constant. When the parameter b on the unit circle this 
constant monotonically increases as arg(6) varies from to ir/2 and diverges when b — > i, 

W« (0) ~ (2vrr ? )- 1 /2 Ekj + Q ^ > b ^ { . (20) 
Other properties of the Faddeev-Volkov model are considered in a separate publication [36]. 
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3 Z-invariant lattices and rhombic tilings 



The above considerations were explicitly carried out for the homogeneous square lattice (shown 
in Fig.l). However, having in mind connections with the discrete conformal transformations 
considered below, it is useful to reformulate the results for more general "Z-invariant" lattices 
[37]. A brief review of this important notion with an account of some subsequent developments 
is given below. 

Following [38] consider a finite set of L directed lines forming a graph XL of the type shown 
in Fig. [3j The lines (in this case six) head generally from the bottom of the graph to the top, 




Figure 3: A set of directed lines going from the bottom to the top of the figure. 



intersecting one another on the way. They can go locally downwards, but there can be no closed 
directed paths in XL. This means that one can always distort XL, without changing its topology, 
so that the lines always head upwards. Now shade alternative faces of XL as in Fig. SJ Form 
another graph 9 by placing a site in each unshaded face, with edges connecting sites belonging 
to faces that touch at a corner. These sites and edges of 9 are represented in Fig. H] by open 
circles and solid lines, respectively. For each intersection of lines in XL, there is an edge of 9 
passing through it, and conversely. The graph XL is the medial graph of 9- Introduce additional 
notation. Let F(S), E(S) and V(S) denote the sets of faces, edges and sites (vertices) of 9 and 
Vi n t (2) the set of interior sites of 9- The latter correspond to interior faces of XL (with a closed 
boundary). 

Now we define a statistical mechanical spin model on 9- With each line £ of XL associate its 
own "rapidity" variable p£. At each site i of 9 place a spin a £ M.. The edges of 9 are either of 
the first type in Fig. [21 or the second, depending on the arrangement the directed rapidity lines 
with respect to the edge. For each edge introduce a "rapidity difference variable" 6 e defined as 

p — q, for a first type edge, 

(21) 

7r — p + q, for a second type edge, 

where p and q are the associated rapidities, arranged as in Fig. [2j Each edge is assigned with 
the Boltzmann weight factor W0 e (a — b), where a, b are the spins at the ends of the edge. 
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Figure 4: The graph S formed by shading alternate faced of Fig. [3]and putting sites on the unshaded 
faces. 



The partition is denned as an integral over all configurations of the interior spins with the 
weight equal to the product of the edge weights over all edges of S , 

z= ( n we^i-*]) n t 22 ) 

(«j)es(S) »eVi„t(S) 

As before, the exterior spins are kept fixed. Taking into account the symmetry (|10p one can 
easily see that the expression (|3j) is just a particular case of ([22]) for the homogeneous square 
lattice of Fig. Q] 




Figure 5: A pictorial representation of the first inversion relation in f|13[) . 

The partition function (|22p possesses remarkable invariance properties [37-39]. It remains 
unchanged by continuously deforming the lines of £ (with their boundary positions kept fixed) 
as long as the graph L remains directed. In particular, no closed directed path are allowed to 
appeao It is easy to see that all such transformations reduce to a combination of the moves 

Actually, these restrictions can be removed if one properly defines "reflected" rapidities for downward going 
lines (see Sect. 3 of [38]), but we will not elaborate this point here here. 
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shown in Fig. [SJ Fig. and Fig. [7J The special role of the star-triangle (|12|) and inversion 




a=b 
o 



p q 



Figure 6: A pictorial representation of the second inversion relation in ([13 





Figure 7: A pictorial representation of the star-triangle relation (j 12|) . 

(|13p relations is that they guarantee the invariance of the partition function under these moves. 
Given that the graphs £ and 9 can undergo rather drastic changes, the above "Z- invariance" 
statement is rather non-trivial. 

The partition function (|22p depends on the exterior spins and the rapidity variables 
pi,P2, ■ ■ ■ ,Pl- Of course, it also depends on the graph £, but only on a relative ordering 
(permutation) of the rapidity lines at the boundaries and not on their arrangement inside the 
graph. Naturally, this graph can be identified with an element of the braid group. Then the 
partition function Z can be regarded as the corresponding representation matrix, acting on the 
continuous boundary spins. There are several celebrated appearances of the braid group in the 
theory of integrable systems, particularly, in connection with the Jones polynomials [39] and 
the Tsuchiya-Kanie monodromy representations in conformal field theory [40]. In Sect. 5 we 
will show that the rapidity graph £ also describes combinatorial properties of integrable circle 
patterns. 

The integral (I22p is well defined in the physical regime, when all the rapidity differences d(ij) 
lie in the interval < Ouj) < 7r. It rapidly converges due to the fast exponential decay (fTTj) of 
the edge weights. For other values of 0(ij) the partition function (I22p is defined with an analytic 
continuation from the physical regime. In what follows we will restrict ourselves to the physical 
regime only (in this case the graph £ cannot contain more than one intersection for the same 
pair of lines). 

Let the graph S has a very large number of edges N. We assume that the number of its 
exterior sites scales as and that the boundary spins are kept finite. Following [37] one can 
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show that the leading asymptotics of the partition function (|22|) at large N has the form 

logZ = - Yl fedge(0 {ij) ) + O(VN) (23) 

where the function f e dge(G) does not depend on specific details of the lattice. It will be the 
same as for the regular square lattice. This result holds for any Z-invariant system with positive 
Boltzmann weights for a large graph XL in a general position. Then, it follows from (I18p . that 
in the Faddeev-Volkov model with our normalization of the weights © the edge free energy 
vanishes identically, 

&) = 0. (24) 

Consider some additional combinatorial and geometric structures associated with the graph 
XL. First, if the unshaded faces in the above definition of 9 are replaced by the shaded ones, 
one obtains another graph 9*, which is dual to 9- Each site of 9* corresponds to a face of 9 
and vice versa. Obviously, both graphs 9 and 9* have the same medial graph XL. Assign the 
difference variables 9 e * to the edges of 9* by the same rule (f2~T]) . Note that there is one-to-one 
correspondence between the edges of 9* and 9- Moreover, if e S E(5) is of the first type then 
the corresponding edge e* G E(S*) is of the second (and vice versa). In other words for the 
corresponding edges 6 e + 9 e * = tt. Let star(i) denote the set edges meeting at the site i. One 
can show that for any interior site of 9 

%) = 2tt, i € V int (S) ■ (25) 

(ij)Sstar(i) 

The similar sum rule holds for the dual graph 9*, 

^2 e {kl) =27r, keV int (S*) (26) 

(kl)Gstar(k) 

There is a dual statistical mechanics model with spins a* attached to the sites i E V(S*)- 
Its partition function is defined as 

n w %j<-°D n w 

(ij)e-B(S*) i<ZVint(9*) 

where the spins at the exterior sites of 9* are fixed. Evidently, it possesses exactly the same 
invariance properties under the deformations of XL, as the partition function (|22p . Using the 
standard arguments [41] based on the duality transformation (|14p one can relate (|27p with (|22p . 

Consider yet another graph XL*, dual to XL. The set of sites of XL* consists of those of 9 and 9*- 
They are shown in Fig. [8] by white and black dots, respectively. The edges of £* always connect 
one white and one black site. The faces of £* corresponds to the vertices of XL. Since the latter 
are of the degree four, all faces of XL* are quadrilateral. The edges of 9 and 9* are diagonals 
of these quadrilateral (see Fig. [9]). Remarkably, the graph XL* admits a rhombic embedding into 
the plane. In other words this graph can be drawn so that all its edges are line segments of the 
same length and, hence, all its faces are rhombi, as shown in Fig. El The corresponding theorem 
has been recently proven in [421. It states that such embedding exists if and only if (a) no two 
lines of XL cross more than oncqj and (b) no line of XL crosses itself or periodic. Note, that in the 
physical regime these conditions are obviously satisfied. 

2 The lines of L are called "train tracks" (!) in [42]. 
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Figure 8: A rhombic embedding of the graph L* 



Assume the edges of the quadrilaterals to be of the unit length and consider them as vectors 
in the complex plane. To precisely specify a rhombic embedding one needs to provide angles 
between these vectors. A rapidity line always crosses opposite (equal) edges of a rhombus. 
Therefore, all edges crossed by same rapidity line p are equal to each other. They are determined 
by one vector £ p , \( p \ = 1. Thus, if the original rapidity graph £ has L lines, there will be only 
L different edge vectors. Choose them as ( Pk = e~ Wk , where pi,p%, ■ ■ ■ ,Pl are the corresponding 
rapidity variables. Each face of £* is crossed by exactly two rapidity lines pk and p£. To this 




Figure 9: Two types of rhombi. 

face associate a rhombus with the edges Q Pk and ( Pe , as shown in Fig. [9J Its diagonals are edges 
of 9 and 9*- The rhombus angles are precisely the "difference variables" 6 e and 9** assigned to 
these edges (this is true for both types of rhombi shown in Fig. [9]). Glue two rhombi together if 
they are adjacent. In the physical regime all their angles will be in the range from to ir. So 
the rhombi will have a positive area and, thus, will not overlap. The sum rules (1251) and (|26H 
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guarantee that the resulting surface is flat with no cusps at the sites of £*. 



4 Quasi-classical limit 



Consider the limit r/ — > oo which corresponds to the quasi-classical limit of the model. The 
parameter rj becomes large when b — > or b — > oo. For definiteness assume b — ► 0, then 
rj ~ (26) . The weight function W$(s) acquires a very narrow bell-shaped form and rapidly 
decays outside a small interval \s\ < b/8. The integral (|22|) can be then calculated with a saddle 
point method. The asymptotic expansion of the Boltzmann weight ([6]) reads 



exp 



1 



2nb 2 



A{e\p)+B{6\p) + 0{b 2 ) \ , 



, 



where 



A(6\p) = A(0\-p) 



log 



1 + e' 



(28) 



(29) 



Note that > 0. It is a smooth even function of p with the minimum value ^4(0|O) = 

reached at p = 0. It is simply related to the Euler dilogarithm function 



A{6\p) =iLi 2 [ -r 
The constant term in (i28|) reads 

Estimating the integral (|22p for small b one gets 



where 



log(l - x) 



9 „,„, \ 1 



zdz 



sm z 



logZ = -^A[p^] + B[^)] + 0(6 2 ) , 



A\p}= Yl A ( e (ij)\pi- Pi) 



and 



i 



log det 



d 2 A[p] 



dpidpj 



&0e.E(S) 



(30) 

(31) 

(32) 
(33) 

(34) 



Here we use rescaled spin variables p = {pi, p 2 , ■ ■ ■}, given by pi = 2-nbui. The symbol p^ 
denotes the stationary point of the action (|33p . defined by the classical equations of motion 



9A[p] 



dpi 



p=p 



(cZ) 



0, i =€ V int (5) 



Explicitly, they read 



n 



(ij)dstar(i) 



Qpj _|_ QPl + W(lj) 
Qpi _|_ Qpj+' e (ij) 



i =€ K„t(S) 



(35) 



(36) 
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Figure 10: A star four edges on the square lattice. The ■(/'-variables are placed in the sites of the 
dual lattice, located in the shaded areas. 

where the product is taken over the star of edges connected to the site i. These are the so-called 
cross ratio equations [22]. They arise in connection with the circle patterns [22], which will be 
considered in the next section. 

The cross ratio equations (j36|) are closely related to the Hirota difference equation [23]. As 
an illustration consider the case regular square lattice S shown in Fig. Q3 A typical "star" of 
four edges is shown in Fig. [TUJ The corresponding equation in (|36p takes the form 

' qPW _|_ gP+i0\ f qPN _|_ qP+\8* \ / qPE _|_ gP+i6>\ /gPS_|_ e P+ifl 



(37) 

gP _|_ qPw+'9 I \ gP _|_ qPn+i9* / \ gP _|_ qPe+i9 / I gP _|_ gps+i^* ' 

where 6 = p — q and 6* = it — 9. Now place a purely imaginary variable ipi, < Imipi < 2ir, 
on every site i of the dual lattice. These sites are shown by black dots located in shaded areas 
in Figs. [10] and [TTJ Connect these new variables to the existing variables pi by the following 
relations. Let i^n and V's be the variables located above and below of a horizontal edge, as in 
Fig. [TTJ Then, we require that 

gPW _|_ ppE + 'lO 

g^iv-^s = M .„ . (38) 

qPE -\- QpW + IV y ' 

Similarly for a vertical edge, 

e fe-^ = e + e ^ e PN- PS = e + e (39) 

gPs 4- gpjv+i» g^B + qv>w+iv v y 

The consistency of these definitions across the lattice is provided by (|37|) . Note that the second 
form of (|39p is identical to f|38|) upon interchanging all ip- and p- variables. Therefore it is natural 
to associate this universal equation with every face of the graph XL* (there are two type of faces, 
as shown in Fig. [9]). 

This the famous Hirota equation [23] , which was the starting point of various considerations 
in [24]. A special feature of the present case is that the p- variables are real p E R, while the 
ip- variables are purely imaginary, |e^| = 1. 
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Figure 11: Faces of the graph S* associated with the equations (|38|) and ([39 



4.1 Star-triangle relation 

Before concluding this section let us discuss some consequences of the Z-invariance. Obviously, 
this property applies to every term of the quasi-classical expansion ([32]) , As an example, consider 
the star-triangle relation (|12p . Introduce new variables 



6>i = 7r - q + r, 2 = 7T — P + q, 9 3 =p — r (40) 

such that 

#1 + #2 + #3 = 2vr. (41) 
Substituting ([28]) into (fT2|) and using the symmetry (fT0|) one obtains 



/ ^Tb ew {-2^ A * [p]+ ®+ [p] } = exp fi AA[pl+13A[/,l+0(h2) } (42) 

where 

A+\p] = A(6i\ Po - pi) + A(6 2 \p - p 2 ) + A{9 3 \p - p 3 ) , (43) 

A A [p] = A(tt - 0i|p 2 - Ps) + A(tt - 2 \p 3 ~ Pi) + A(ir - 9 3 \p x - p 2 ) ■ (44) 

Expressions for and "B/\ are defined in a similar way, with the function A(9\p) replaced by 
B(9\p), given by pip. The star-triangle relation (|42p implies two non-trivial identities valid for 
arbitrary values of p\ , p 2 , p 3 ■ The first one is 

A+[Pq cl \pi,P2,P3] = A A [p 1 ,p 2 ,p 3 ], (45) 

where p^ is the stationary point of the integral in (|42p 



(d) , / r 2 r 3 sin 6>i + n r 3 sin 6> 2 + H r 2 sin 6> 3 \ 

Po -log — — — — — , n = e ft . (46) 

\ n sm t^i + r 2 sin # 2 + r 3 sin 9 3 J 

With an account of (|30p it is a "twelve-term" dilogarithm identity, which can be easily verified 
by using the definition of the dilogarithm. A geometric meaning of this identity in terms of 
volumes of certain polyhedra in the Lobachevskii 3-space is explained in the Appendix A. The 
second identity 

p 0=p ^> " L ' 'Po=P 



2 dp 



2 



,.„ - £*[pL_„W-£aM (47) 
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is an algebraic identity for r\ , r-i , r% , which follows from elementary (but lengthy) calculations 
with the use of the explicit expression (|46p and the relation 



Note that exactly the same quasi-classical star-triangle relation (|42p arises in a related, but 
different, context of the chiral Potts model [43]. 

5 Circle patterns and discrete conformal transformations 

The Faddeev-Volkov model on a Z-invariant lattice has a remarkable connection with the circle 
patterns related to the discrete Riemann mapping theorem [12-15]. The classical Riemann 
mapping theorem states that for any simply connected open domain in the complex plane 
(which is not the whole plane) there exists a biholomorphic mapping from that domain to the 
open unit disk. The fact that the mapping is biholomorphic implies that it is a conformal 
transformation. Likewise, the discrete Riemann mapping theorem involves discrete analogs of 
conformal transformations. Below we give a short introduction into this subjectH. 

In the continuous case the conformal transformations possess two key properties: (i) they 
preserve the angles and (ii) uniformly rescale all infinitesimal lengths in a vicinity of every 
point, where the scale depends on the point. Now one needs to reformulate these properties 
for a discrete case, when the continuous plane is replaced, for example, by a polygonal lattice. 
For the first property one could, naively, suggest that the angles between the edges meeting at 
the same site remain unchanged. We will call this the rule (i). The second property is more 
complicated, since there is no a lattice analog of the infinitesimal length. Nevertheless, as a first 
attempt it is reasonable to try the simplest possibility and suggest that the lengths of the edges 
attached to the same site are scaled equally (though the scale can vary from site to site). We 
will call this the rule (ii). 

Now consider transformations of the ordinary square lattice, implementing both rules at 
every site. The rule (i) requires that the 90° angles between the four edges at any site remain 
unchanged. The rule (ii) requires that these edges also remain equal in length for every site 
(though, in principle, the lengths could change from site to site). A very little inspection shows 
these requirements only allow a uniform dilatation of the whole lattice (of course, the same 
would be true for any polygonal lattice, not just the square one). Obviously, the rules are 
too restrictive and do not allow any interesting transformations of the lattice; they need to be 
relaxed. 

A fruitful idea is to work with a bipartite lattice and apply only one rule on each sublat- 
tice. Consider, for instance, the square lattice shown in Fig. [T2l Its sites are divided into two 
sublattices shown with "white" and "black" dots. Denote them as S and S* respectively (to 
unload the picture the edges of S, connecting the white sites, are not shown). There are circles 
centered at every white site and passing through four neighboring black sites. Thus every face 
of black lattice, S*, is inscribed in a circle. Two circles intersect at two distinct points (not 
just touching) only if the corresponding faces are adjacent. In the case of Fig. [12] they always 
intersect at the right angle, which is the angle between the edges normal to the circles. 

Let us now distort our lattice implementing the rule (i) for all black and rule (ii) for all white 
sites. A particular example of such transformation is shown in Fig. 1131 As seen from the picture, 

3 We are indebted to A. Bobenko and B. Springborn for many illuminating explanations, some of which are 
used below. 




(48) 
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Figure 12: The regular square lattice and an associated circle pattern. 



the circle centered at white sites still pass through all neighboring black sites (rule (h)), they still 
intersect at the 90° angles (rule (i)), but the radii of these circles have changed! It is important 
to keep in mind that, by construction, this transformation preserves the combinatorics of the 
original lattice. In particular, there is a one-to-one correspondence between the circles on the 
two pictures, which preserves their adjacency and intersection angles. 

Before formulating equations determining the radii of the circles let us make a few remarks. 
Denote by z = {z\, z%, . . . , zn} positions of the circle centers in the complex plane on the original 
square lattice, where N is the total number of circles. In the case of Fig. [12] this number is 49. 
The circles there cover a square-like region, bounded by a dashed line. Let {z[, z' 2 ,..., z' N } be 
the centers of the corresponding circles in Fig. [T3l The transformed circle pattern covers a disk- 
like region, bounded by a large circle, also shown by the dashed lin43- Now define a function 
/ such that z' k = /(#fc), for all Zk G z. This is a discrete approximation to the continuous 
Christophel-Schwartz map. When N tends to infinity, f{z) uniformly converges to the latter 
with an accuracy 0(N~ l ) [14,15,44]. 

The patterns considered above enjoy the combinatorics the square lattice. All circles there 
intersect at the 90° angles. In a general case, one can consider circle patterns with the combina- 
torics of an arbitrary planar graph S, where the edges are assigned with arbitrary intersection 
angles {6* e }, ee£'(S),0<6' e <7r. The centers of the circles correspond to the sites of 9 (white 
sites), such that every face of the dual graph S* is inscribed in a circle. Two circles connected 
by an edge e £ -E(S) intersect at the angle 9 e assigned to this edge (Fig. [T6j) . Their intersection 

4 The pattern of Fig. [13] is uniquely fixed (up to a uniform rotation) by the four-fold axial symmetry and a 
requirement that all exterior circles there touch one additional circle. 
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Figure 13: A circle pattern with the combinatorics of the square lattice. 



points are sites of S* (black sites). Examples of two different circle patterns associated with 
same graph S as in Fig. [Hand having the same intersection angles are shown in Figs. 1151 and 1141 
Thus, to define the intersection properties of a circle pattern (its combinatorics and intersection 
angles) one needs to specify a set {S, {# e }}- 

The intersection properties are called integrable [45] if and only if the corresponding cir- 
cle pattern admits an isoradial realization (where all its circles are the same). Here we will 
only consider this case. The integrable circle patterns were extensively studied in connection 
with various approximation problems for the Riemann mappings which are typically based on 
a regular lattice combinatorics [15]. They also connected with integrable discrete non-linear 
equations [45]. 

The space of all "integrable sets" {S,{#e}} was completely described in [45]. Essentially the 
same arguments were already presented in Section 3 throughout the discussion of the Z-invariant 
systems. Remind that our starting point there was a rapidity graph XL (Fig. [3]) with a set of 
rapidity variables, one per each line of XL. Using these data we constructed a rhombic "tiling" 
of the dual graph XL* (Fig. [8]). From this tiling one can immediately construct an isoradial 
circle pattern by drawing circles of the same radius centered at the white and passing though 
the neighboring black sites (Fig. [T4"|) . Conversely, any isoradial circle pattern defines a certain 
rhombic tiling from which one can easily reconstruct [42] the corresponding rapidity grapfjl 

5 Note, that in general this graph can have a multiply connected boundary. 



16 



Figure 14: An isoradial circle pattern obtained from the rhombic tiling of Fig. [8] The symbol <E> 
denotes the cone angle at an exterior site. 



together with all rapidity values (up to an insignificant overall shift). Thus the intersection 
properties of any integrable circle pattern can be completely described in terms of some rapidity 
graph, arising in the Z- invariant models (and vice versa). 

Let us now derive the equations determining the radii of the circles. We assume the same 
notations as in Section 3. We consider the circle pattern with the combinatorics of 9 and identify 
the rapidity difference variables &Uj), (ij) G E($) with the circle intersection angles. Remind 
that these variables obey the sum rules (|25p and (|26p . Every pair of intersecting circles with the 
centers i and j, i,j £ 9, and the radii rj and rj is associated with a kite-shaped quadrilateral 
shown in Fig. [T6l The angles ipjk and ipjk bisected by the edge (jk) are given by (in the notations 
defined on the figure) 

<Pjk = ~ log — z^- , (49) 

% r j + r k e 

and ipkj is obtained by permuting rj and r^. Consider a circle flower consisting of one central 
circle and a number adjacent circles (petals), like the one shaded in Fig. [P31 Obviously the kite 
angles at the center of the flower add up to 2n 

^2 na) = 2n > ie v ^t(s) ■ (so) 

(ij)(zstar(i) 



Similarly, a sum of the kite angles at any interior black site (these are 9* angles as in Fig (16 
should be equal to 2ir as well, but this condition is automatically satisfied in virtue of (|26p 
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Figure 15: An integrable circle pattern with combinatorics of the graph 3 from Fig. |U The sum of 
the kite angles at the center of the circle flower is equal to 2ir. 




Figure 16: Two intersecting circles with the centers at the sites j and k. The angle 9jk is assigned 
to the edge connecting these sites. 



Introduce logarithmic radii pi = log rj of the circles. Substituting (fl9|) into ([50]) and taking 
into account the sum rule (|25H one precisely obtains the equation ([36]) which arose in the quasi- 
classical limit to the Faddeev-Volkov model. 

The quasi-classical action A[p] in (|33p correspond to the fixed boundary conditions, where the 
radii of all boundary circles are fixed. They can be chosen arbitrarily. Bobenko and Springborn 
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[27] considered free boundary conditions. Their action 



A B s[p] =A[p\- ^2 Pi®i + const 

ieV ext (9) 



(51) 



involve cone angles <3?j at the centers of exterior circles (see Fig. [T4")) . Their values are arbitrary 
modulo a single constraints 



where Ni n t is the number of interior sites of 9- In the variation principle for the action (|5ip one 
considers the angles <3?j as fixed and the external radii as independent variables. 

6 Conclusion 

In this paper we displayed some remarkable connections of the theory of integrable quantum 
systems to discrete geometry. We have shown that the Faddeev-Volkov model provides a quanti- 
zation of the circle patterns and discrete conformal transformations associated with the discrete 
Riemann mapping theorem. 

The partition function of the model in the thermodynamic limit is calculated here via the 
inversion relation method. It would be interesting to verify this result by other means (par- 
ticularly, by numerical calculation) since the inversion relation method is based on unproved 
analyticity assumptions of the "minimal solution" ([8]) . At the moment not much else is known 
about the Faddeev-Volkov model. In particular, the problem of the diagonalization of the trans- 
fer matrix is not yet solved, though some important advances in the related problems in the 
continuum sinh-Gordon were made in [46,47]. 

Finally mention the intriguing relation of the Faddeev-Volkov model with the hyperbolic 
geometry discussed in the Appendix A. The solution ([6]) is constructed from the non-compact 
quantum dilogarithm @. Note, that there exist another solution of Yang-Baxter equation 
connected to this dilogarithm [48]. The latter is related to the link invariants [49] which are 
also connected with hyperbolic geometry, particularly, to the volumes of hyperbolic 3-manifolds 
[50,51]. It would be interesting to understand these connections further. 
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Appendix A. Star-triangle relation and hyperbolic geometry 

An excellent introduction into the volume calculation in hyperbolic geometry can be found 
in [52,53]. As shown in [54,55] the action (j33H has a geometric interpretation as the volume 
of a certain polyhedron in the Lobachevskii 3-space. Consider the Poincare half-space model 
{x,y,z E M|z > 0} of the hyperbolic 3-space with the metric ds 2 = (dx 2 + dy 2 + dz 2 )/z 2 . In 
this model hyperbolic planes are represented by hemispheres and half-planes which intersect the 
infinite boundary orthogonally in circles and straight lines. Take a circle pattern (for instance, 




(52) 



i£V ext (S) 



e€E(S) 
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the one shown in Fig. I14|) and imagine it as lying in the infinite boundary (the xy-plane). If we 
erect hemispheres over the circles and orthogonal half-planes over the prolonged boundary edges 
(these edges are shown by heavy solid lines in Fig. I14|) . we obtain a set of hyperbolic planes which 
bound a polyhedron. Most of its vertices are at infinity; these include the circle intersection 
points (in the xy-plane) and one extra point at the infinite boundary — the intersection point 
of the planes raised from the boundary edges. Obviously, this polyhedron can be glued from 
similar shaped hyperbolic hexahedrons, shown in Fig. [T71 which are associated with every pair of 
intersecting circles. Every such hexahedron can be split into two tetrahedra with three vertices 
at infinity. 




Figure 17: A hexahedron in the hyperbolic 3-space bounded by four infinite vertical planes and 
two hemispheres. 

In what follows we will refer to Fig. 1161 which corresponds to the projection of the hexahedron 
to the ideal, but simplify the notations there taking j = 1, k = 2 and denoting (pj^, ip^j, 9* k 
by 0i, 02, O, 0*, respectively. Note that cfii + <p2 = 0. From ff^9|) one has 

pPl _L P p2+i0 

= +e . (A.l) 
where pi = logrj. Introduce the Milnor's Lobachevskii function 

Jl(x) = - [ log|2sin£|d£ . (A.2) 



The volume of the hexahedron in Fig. [T7] is given by [54] H 

V hexahedron (<Pl,<P 2 ) = Jl(y) + Jl(y) = 2Jl(-) + ^A(9\pi-p2) + 2(<h-<h)(fil-P2), (A.3) 



'Note, that our variables 4>i,cf>2,0 correspond to 2(j>\,2(j>2,9* in [54]. 
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where the function A(Q\p) is defined in (|29|) . Then, the quasi-classical action yi[p( c ^] in (|32|) 
(remind that it corresponds to the fixed (Dirichlet) boundary conditions) is given by 

A[pW] = 2V(P)-2V(P )- E ( A - 4 ) 

i&Vext(S) (ij)Gstar(i) 

where V(P) is the volume of the hyperbolic polyhedron P whose projection to the ideal is a 
circle pattern with a given set of intersection angles 0« and boundary radii rj = e Pi . Further, 
V(Pq) denotes the value of this volume for the isoradial case (when all circles have the same 
radius) 

V(P )= £ Jl( 9 f). (A.5) 

(ij)eE(S) 

Thus, up to boundary terms the action (|A.4p coincides with an "excess volume" of P with respect 
to Pq. The result (|23p and (|24|) implies that if the polyhedron P is glued from a very large 
number of hexahedrons, N, this excess volume grows only as \/N. This is not very surprising: 
an arbitrary integrable circle pattern is significantly different from its isoradial counterpart only 
at the boundary (indeed, compare Fig. [TJ] and Fig. [15]) , 




Figure 18: Circle patterns corresponding the star (left) and triangular sides of the Yang-Baxter 
equation ([42]) . 

The star-triangle equation has the following interpretation. The volume of the polyhedron 
erected over circle pattern corresponding to the star side of (142p . shown in Fig. [TU is equal to 

V+ = 2Ji(^) + 2Jl(^) + 2Jl(^) +^A ir [p { cl \ Pl ,p2,P3]+ boundary term (A.6) 

where the action and p^ are defined in (|4"3"|) and (|4"o]) . Similarly, the volume of the polyhe- 
dron corresponding to the triangle side of the star-triangle relation in Fig. [18] is given by 

V A = 2Jl(^^ i ) + 2Jl(^-^) + 2Jl(^-^) + ^ A bi,P2,p 3 ] + boundaryterm (A.7) 

where A a is defined in (|44p . Clearly, the difference between these two volumes is the volume of 
an ideal tetrahedron with its vertices located at the circle intersection points in Fig. [TBJ Indeed, 
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taking into account (|45|) and the fact that the boundary terms in (|A.6|) and (|A.7p coincide, one 
obtains 

tetrahedron = V* ~ V A = 2 Jl(|) + 2 Jl(|) + Jl(|) - 2 Jl(^i) - 2 Jl£^) - Jl(^^) . 

(A.8) 

Using the identity 

2JI(^)-2JI(^) = JI(0). (A.9) 

one obtains 

tetrahedron = Jl(fll) + Jl(0 2 ) + Jl(0 3 ) , (A.10) 

which is precisely the Milnor's formula [52] for the volume of the ideal tetrahedron. We would like 
to stress that whereas the volumes V+ and V/\ depend on the radii of the circles, their difference 
does not. It depends only on the angles 9\, 62 and #3. In (f4"5j) this difference is absorbed into 
the normalization of the edge weights. Correspondingly, the RHS of (IA.8j) contains six terms, 
one for each edge of the tetrahedron. Clearly, the action (|A,4|) is invariant under any star- 
triangle transformation since volumes of ideal tetrahedrons will cancel out from the difference 
V(P)-V(P ). 



Appendix B. Properties of the functions cp(z) and <fr(z). 

Below we use the following notations 

»?= 5(6 + 0, q = e- b \ q = e-™-\ q = \e*p ^-^ , (B.l) 

00 

(x,q)oo = 11(1 -g*x) . (B.2) 

The function ip(z). This function is defined by the integral (|4|). It has the following properties. 

(a) Simple poles and zeros 

poles of tp(z) = {ir/ + \mb + \nb~ l , m, n G Z>o} , 
zeros of 9? (z) = { — ( 177 + \mb + \nb~ l ) , m,nSZ>o} . 

(b) Functional relations 

v(zM-z) = e'^-a-V)/- , = (l + e*™* 1 ) . (B.4) 

(c) Asymptotics 

<p(z) ~ 1, Re(z) -00; p(z) ~ e i ^ 2 -™ (1 - 2??2)/6 Re(z) +00 . (B.5) 
where Im(z) is kept finite. 

(d) Product representation 

(—np 27rzb - a 2 ) 

?(*) = t 2^ ~t 1 Imfo2 > ( B - 6 ) 
(-qe 2nzb ;q 2 )oo 
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(e) Pentagon relation. The function ip(z) satisfy the following operator identity [3] 

p(PMX) = ^(X)^(P + X)^(P) , [P,X] = -L . (B.7) 

Z7TI 

where [ , ] denotes the commutator. It can be re-written in the matrix form [56, 57] 

g(g + u ) e 2ir\wx dx = e i7T(l+V)/12 -2wiw(v+\ V ) ip(u - V - \1])<p(w + \ri) ^ 

(p(x + v) ip(u — v + w — \rj) 

where the wedge of poles of (p(x + u) must lie in the upper half-plane, the wedge of zeros of 
ip(x + v) must lie in the down half-plane, and the integrand must decay when x — ► dboo. 

The function <&(z). This function is defined by the integral ©. It has the following properties. 

(i) Simple poles and zeros 

poles of &(z) = {2iry + \mb + \nb~ l , to, n £ Z>o , to + re — \m — n| = mod 4} , 
zeros of $(2) = { — (2ir/ + \mb + ireft -1 ) , to, re 6 Z>o , m + n — |m — n| = mod 4} . 

(B.9) 

(ii) Functional relations 

*(*)$(-*) = e^ 2 / 2 "^ 1 ^ 2 )/ 12 , *(z + ir/)$(z - ire) = ^(z) . (B.10) 

(iii) Asymptotics 

*(z)~l, Z^-OO; ~ e i7r ^/2-i 7 r(l-8'7 2 )/12 ; z _ +QO _ ( BU ) 

where Im(z) is kept finite. 

(iv) Product representation 

f 2 2m6. 4\ / 77^-^^/(217) . — 2n 

^/ \ w e ) 1 ho \ ae i y Joo T ,2 . n /D 10 \ 

*W- M v.^ it 4). B .-vw>,^ ' Imfc>0 - (R12> 
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